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ABSTRACT 

We present new radio observations of the supernova SN 1979C made with the VLA at 20, 6, 3.6, and 
2 cm from 1991 July to 1998 October, which extend our previously published observations (Weiler et 
al. 1986, 1991), beginning 8 days after optical maximum in 1979 April and continuing through 1990 
December. We find that the radio emission from SN 1979C has stopped declining in flux density in the 
manner described by Weiler et al. (1992), and has apparently entered a new stage of evolution. The 
observed "flattening," or possible brightening, of the radio light curves for SN 1979C is interpreted as 
due to the SN shock wave entering a denser region of material near the progenitor star and may be 
indicative of complex structure in the circumstellar medium established by the stellar wind from the red 
supergiant (RSG) progenitor. 

Subject headings: stars — stars: evolution — supernovae: individual (SN 1979C) 



1. INTRODUCTION 

The study of supernovae (SNe) which are significant 
sources of radio emission, known as "radio supernovae" 
(RSNe), provides unique information on the properties of 
the progenitor stellar systems and their immediate circum- 
stellar environments. In particular, changes in the density 
of the presupernova stellar wind established circumstellar 
material (CSM) alter the intensity of the radio emission 
and allow us to probe the mass-loss history of the super- 
nova's progenitor, structures in the CSM, and the nature 
and evolution of the SN progenitor. 

Significant deviation of the radio emission of SNe from 
standard models has previously been noted and inter- 
preted as due to a complex CSM density structure. 
SN 1979C (Weiler et al. 1986, 1991) has shown a quasi- 
periodic variation in its radio emission which may be due 
to modulation of the CSM density by a binary companion 
(Weiler et al. 1992; Schwartz & Pringle 1996). SN 1987A, 
after an initial, faint radio outburst and rapid decline un- 
til 1990 June, is now increasing in radio flux density due 
to the SN shock starting to impinge on the inner edges of 
the well known, much higher density central ring (Turtle 



et al. 1990; Staveley-Smith et al. 1992, 1993, 1995; Ball 
et al. 1995; Gaensler et al. 1997). SN 1980K (Montes et 
al. 1998) at an age of ~ 10 years has experienced a sharp 
drop in its radio emission far beyond that expected from 
models of its previous evolution. More recently, SN 1988Z 
has also shown a sharp drop in its radio emission similar 
to that of SN 1980K (Lacey et al. 1999). 

Recent observations of SN 1979C imply that the shock 
from SN 1979C has entered a new, higher density CSM 
structure different from the modulated decline previously 
reported (Weiler et al. 1992). The radio emission has ap- 
parently stopped declining and has been constant, or per- 
haps increasing, for the past eight years. 

We have monitored SN 1979C [RA(2000.0) = 
12h22'"58!67 ± ffiOl; DEC(2000.0) = -hl5°47'51'.'8 ± 0'.'2] 
in NGC 4321 (M 100) since 1979 April 27, and the re- 
sults through December 1990 are available in Weiler et 
al. (1986, 1991, 1992). Here we present new Very Large 
Array^ (VLA) radio measurements between 1991 July 28 
and 1998 October 19. 



^Affiliated with the Astrophysics Division, Space Science Department of ESA. 
^Now at the Massachusetts Institute of Technology. 

^The VLA is operated by the National Radio Astronomy Observatory of the Associated Universities, Inc., under a cooperative agreement 
with the National Science Foundation. 
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2. OBSERVATIONS 

Routine monitoring observations of SN 1979C have been 
carried out with the VLA at 20 cm (1.465 GHz) and 6 
cm (4.885 GHz) a few times a year since 1991. On 1996 
December 21 & 1998 February 10, observations were ad- 
ditionaUy made at 3.6 cm (8.435 GHz) and 2 cm (14.965 
GHz). The techniques of observation, editing, cahbration, 
and error estimation are described in previous pubhca- 
tions on the radio emission from SNe (see, e.g., Weiler 
et al. 1986), with particular reference to SN 1979C in 
Weiler et al. (1991, 1992). The "primary" calibrator re- 
mains 3C 286 and is assumed to be constant in time 
with flux densities 14.45, 7.42, 5.20, and 3.45 Jy, at 20, 
6, 3.6, and 2 cm, respectively. The "secondary" calibra- 
tor, 1252-1-119, was used as the phase (position) calibra- 
tor with a defined position of RA(1950)=12'^52"^07!724, 
DEC(1950)= -hll°57'20'.'82, and, after calibration by 
3C 286, as the actual flux density reference for the ob- 
servations of SN 1979C. As expected for a "secondary" 
calibrator^, the flux density of 12524-119 has been varying 
over the years, as can be seen in both Table and Fi gure 

The flux density measurement error is a combination 
of the rms map error, which measures the contribution of 
small unresolved fluctuations in the background emission 
and random map fluctuations due to receiver noise, and 
a basic fractional error e included to account for the nor- 
mal inaccuracy of VLA flux density calibration (see, for 
example, Weiler et al. 1986) and possible deviations of 
the primary calibrator from an absolute flux density scale. 
The final errors (cr/) as listed in Table || are taken as 

aj EE [eSof + ct'o (1) 

where 5*0 is the measured fiux density, do is the map rms 
for each observation, and e = 0.05 for 20, 6, and 3.6 cm, 
and e = 0.075 for 2 cm. 

3. RESULTS 

Table ^ shows the new fiux density measurements from 
1991 July 28 through 1998 October 19, while Figure | 
shows a plot of all the available 20 and 6 cm data for 
SN 1979C from 1979 Aprfl 27 through 1998 October 19, 
along with the best-fit model through 1990 December (day 
- 4,300; Weiler et al. 1992). Figure | shows the evolution 
of the spectral index between 20 and 6 cm, af*, for all 
the available data through 1998 October 19, and Figure^ 
shows the spectrum at four frequencies for the 1996 De- 
cember 21, 1998 February 10, and February 13 observa- 
tions. 

The solid line in Figure ^ is calculated from the best- 
fit model to the data through 1990 December and is not 
adjusted for the change in evolution since that time. The 
solid line in Figure ^ is calculated for the best-fit spectrum 
using only data from the 1996 December 21, 1998 February 
10, and February 13 observations, and may be represented 
as Si, cx t^" with a = —0.63 ± 0.03. This spectral index 
value is clearly fiatter than both the best-fit spectral index 
after day 4,300, 004300 ^ —0.70, and the pre-1991 best-fit 
spectral index of —0.75 (see Table 3). Even though for the 

''Secondary calibrators are chosen to be compact and unresolved by 
and variable so that their flux density must be recalibrated from the 



new measurements since 1990 December ranges from 
about — 1 ^ <; —0.6 (Figure ||), the spectral index 
retains a completely optically thin, nonthermal character, 
at least over the wavelength range from 20 to 2 cm (Fig- 
ure ^). It is apparent that, while the evolution of the 
radio emission (and presumably the density structure of 
the CSM) has changed significantly since 1990 December, 
the best-fit spectral index, since becoming optically thin, 
has remained fairly constant with at<43oo ~ —0.75 (be- 
fore day 4,300; Figure ||) and Q;t>4300 ^ —0.70 (after day 
4,300), although short-term deviations such as that seen 
in Figure ^ (a = —0.63) appear to be real. 

The fluctuations seen in Figure || in the 6 cm data (rel- 
ative to the fairly constant 20 cm data) after day 4,300 
also appear to be real. This data was reduced indepen- 
dently by three (and in some cases four) of the authors, 
all with similar results. The fluctuations are not correlated 
with fluctuations in the secondary calibrator 1252-1-119 or 
VLA array conflguration. Higher frequency observations 
which would assist in the interpretation of this data are, 
unfortunately, available on only two dates (1996 December 
21 & 1998 February 13) and have no long-term monitoring 
values with which to compare. 



4. PARAMETERIZED MODEL 

Previous work on RSNe (see, e.g., Weiler et al. 1986, 
1991) has shown that the radio emission can be reasonably 
well described in its gross properties by a parameterized 
model of the form 

where 




where Ki and K2 correspond formally to the flux density 
(in mjy) and uniform absorption, respectively, at 5 GHz 
one day after the explosion date, to; a is the nonthermal 
spectral index of the synchrotron emission; and (3 is the de- 
cline rate of the radio emission after maximum. The term 
e~'^ describes the attenuation of a local, external medium 
that uniformly covers the emitting source ("uniform ex- 
ternal absorption" ) and is assumed to be purely thermal, 
ionized hydrogen with frequency dependence From 
the Chevalier (1982a,b) model for radio emission from 
SNe, this CSM is assumed to have radial density depen- 
dence p oc and to have been established by a constant 
mass-loss rate, M, constant speed wind, w, from a red su- 
pergiant progenitor. The parameter S describes the time 
dependence of the optical depth for this local, uniform 
medium, and S = a — (3 — 3 is specified in the Chevalier 
model (Chevalier 1984). For an undecelerated SN shock, 
S = —3 is appropriate (Chevalier 1982a). 

Chevalier (1982a) also determined the dependence of the 
radio luminosity on the mass-loss rate (M) and progeni- 
tor wind velocity (w) or, equivalently, to the average CSM 

longest VLA baselines. Such compact objects arc usually extragalactic 
nary calibrators for each observing session. 
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-7+12m)/4 



(4) 



where 7 = —2a + 1 is the power law of the relativis- 
tic electron energy distribution and m = —S/3 describes 
the time dependence of the self-similar evolution of the 
shock radius, R oat™- . In this model, the magnetic energy 
density and the relativistic energy density both scale as 
the total post-shock energy density (oc PcSM^jhock where 
Kihock = dR/dt), and both the magnetic field amplification 
efficiency and the particle acceleration efficiency remain 
constant as the SN evolves. 

Weiler et al. (1992) noted periodic features in the radio 
light curves of SN 1979C. Since the spectral index was not 
affected by the observed flux density oscillations, they con- 
cluded that the variations were due to emission efficiency 
changes caused by modulations in the CSM density struc- 
ture, rather than being due to optical depth effects. They, 
therefore, introduced a modification to Equation^ by mul- 
tiplying the emission term Ki by the sinusoidally-varying 
modulation 



1 + A sin 
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1 day 
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-(a+2i5+3)/2 



(5) 



density (pcsM oc M/w), as parameters from the actual data. The ensemble of param- 

eter fits is then used to estimate errors for the parameters 
by examining number distributions for the parameter in 
uestion. The errors in the fitting parameters in Table 
correspond to the values with 15.85% and 84.15% (i.e., 
±lcr for a Gaussian distribution), respectively, of the cu- 
mulative distribution for each parameter. 

5. DISCUSSION 

5.1. Interpretation 

As can be seen in Figure |^, the new measurements for 
SN 1979C since 1991 July 28 do not fit the model which 
described the data relatively well through 1990 December. 
Since that time, the fiux densities at all frequencies ap- 
pear to have stopped declining and perhaps even started 
increasing. Accepting our previous interpretation that a 
constant spectral index indicates an unchanged emission 
mechanism with constant efficiencies, this flattening im- 
plies that the shock wave of the SN is now interacting 
with a higher density structure in the CSM. 

Under this scenario, the flattening of the radio light 
curve may imply eithe r the presence of denser circum- 
stellar material (§ 5.2.1 ), or perhaps a geometrica l solu - 
tion, such as the presence of a flared disc or ring (§ 5.2.2 ). 
Such dense shells or rings, similar to those we propose, are 
present around early type stars, as indicated by observa- 
tions of luminous blue variable (LBV) stars (see, e.g., Nota 
et al. 1995), the recent detections of rings around B type 
supergiants (see, e.g., Brandner et al. 1997), and even by 
SN 1987A with its prominent circumstellar ring. 

Another example of a SN surrounded by a circumstel- 
lar shell denser than a normal RSG wind is provided by 
SN 1978K, for which Montes et al. (1997) have found ev- 
idence of a constant free-free absorption that produces a 
characteristic low-frequency curvature in the typical syn- 
chrotron radio spectrum of the SN. The presence of a dis- 
crete shell of ionized gas surrounding SN 1978K has been 
confirmed by Chu et al. (1999) through spectroscopic mea- 
surements showing narrow Ha and [N II] emission lines 
with widths and intensities characteristic of the environ- 
ment around LBV stars. 

Additionally, there is evidence from optical observations 
of SN 1979C that the Ha and [O I] emission lines are 
also not decreasing, but have remained fiat or perhaps in- 
creased slightly (Fesen et al. 1999). The optical emission 
arises from a different emission mechanism than the ra- 
dio emission and at a different region of the SN shock, 
providing a separate line of evidence that SN 1979C has 
encountered a higher density region of the circumstellar 
medium. 

Furthermore, recent results from simulations of young 
SNRs interacting with higher density clouds of varying 
sizes in an otherwise uniform CSM (Jun & Jones 1999 
and references therein) have shown that radio emission in- 
creases due to the increased ambient density and turbulent 
magnetic field amplification after a SNR/cloud interaction 
begins. The simulations also showed changes in the effi- 
ciencies over time arising from the SNR/cloud interaction, 

^While it is true that the optical depth is due to the integral along the line of sight to the radio emitting region, most of the absorption 
occurs close to the shock front where the density is the greatest. The fractional errors in r from using our equations (3) and (6) instead of the 



Since a radial density modulation of the CSM will also 
affect the optical depth t (t oc p^SM °^ (M/w)^), Equa- 
tion ^ was modified by multiplying the absorption term 
K2 by 



l-f A sin 
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In Equations 



1 day 
and ^ the expression 
to 
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A sin 
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1 day 



C 
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(7) 



represents the deviation of M/w from a constant pre-SN 
wind mass-loss rate^. The parameters A,B, and C de- 
fine the sinusoidal variation of pcsM (oc M/w), where A 
represents the fractional amplitude of the density modula- 
tion, B represents its frequency (in cycles day~^), and C 
represents its phase lag (in radians). 

The best-fit model parameters to the pre-1991 [t — tp < 
4,300 days) data are listed in Table ||. This fit is slightly 
different from that of Weiler et al. (1992), due to improved 
fitting software, but agrees with that work to within the 
uncertainties. The errors in these new fitting parameters 
were estimated using a bootstrap procedure (Press et al. 
1992). Bootstrap procedures use the actual data sets to 
generate thousands of synthetic data sets that have the 
same number of data points, but some fraction of the data 
is replaced by duplicated original points. The fitting pa- 
rameters are then estimated for these synthetic data sets 
using the same algorithms that are used to determine the 



actual integrated expression for r vary roughly sinusoidally with a period 
T <g 1, thus the fractional errors in e~'^ are much smaller than 2A. 



1/B, and a magnitude ~ 2A. At the epoch we are concerned with. 
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and how, in such interactions, the efficiencies, magnetic 
fields, and densities do not scale in the simple manner de- 
scribed by the Chevalier (f982a,b) model. 

However, it must be noted that a change in the ob- 
served radio emission could be caused by a change in the 
efficiency of particle acceleration without significant CSM 
density change. Since the efficiency of particle acceler- 
ation is not well understood, an efficiency change may 
or may not cause a noticeable change in the spectral in- 
dex. If we take the 77 to be the ratio between the num- 
ber of electrons that are accelerated to the total num- 
ber of electrons, then 77 is linearly proportional to the 
measured flux density. While we may not know 77, we 
do know that the ratio of the measured flux density to 
that expected from the Chevalier model. If we assume 
that only rj changes, while the densities and magnetic 
fields behave as expected from the Chevalier model, then 
'^Tioo days /'^4300 days ^ 'S'mcasurcd/5'inodoi ^ 1-7. This 
would require that the acceleration efficiency had remained 
roughly constant until ~ 4,300 days after explosion, then 
continually increase, being some ^ 70% greater by day 
~ 7,100 than at day ~ 4,300. Nevertheless, given the op- 
tical evidence of an increase in the circumstellar density, 
which does not rely on particle acceleration, we interpret 
the deviation of the radio emission of SN 1979C from the 
previously assumed model as due to local density enhance- 
ment with the particle acceleration efficiency relatively un- 
changed. 

5.2. Nature of the Necessary Density Enhancement 

5.2.1. CSM Density Increase 

A denser CSM may either be an isotropic shell or be in 
the form of condensations which have a small combined 
covering factor cj) = ^/A-k < 1, relative to the expand- 
ing shock front which is overtaking them. Since the flux 
density increase at t ~ 7,100 days is about a factor of 
1.7 greater than that expected from the previous best- 
fit model (Weiler et al. 1991), and since the radio lumi- 
nosity has a power-law dependence on the CSM density 

(-^radio 

°^ PcsM' Chevalier 1982a), the implied density en- 
hancement is a factor of ^ 1.340^"'^^ greater than the 
expected density for a p oc r^^ constant mass-loss rate, 
constant velocity wind established CSM. Moreover, the 
increase relative to the best-fit model behavior is gradual 
with time, indicating a relatively smooth density distri- 
bution, which decreases with radius more slowly than the 
expected r~^ behavior. To produce an increase of density 
by a factor ~ 1.34 between 4,300 and 7,100 days a new 
radial density profile of approximately pcsM cjc r~^'^ is 
required if </) = 1, i.e., for a imiform density distribution. 

If the density enhancement is the result of discrete fea- 
tures, whether a shell, condensations, or a ring, it is likely 
to fill a relatively limited fraction of the space around the 
SN progenitor, and to be present over a limited interval of 
radii (AR/R < 1). Thus, this increase of radio emission 
could be a transient phase that disappears in a relatively 
short time. 

The argument for a limited duration to the increased 
radio emission, at least for a spherical geometry, is based 
on the plausible initial mass of the SN's RSG progeni- 
tor. The best-fit parameters for the radio emission at 
t < 4,300 days (Table 3) imply a pre-SN mass-loss rate 



of M ~ 1.6 X 10~^ Mq yr"\ if one adopts a pre-SN 
stellar wind velocity of 10 kms~^, a shock velocity of 
9,250 km s~^ (Bartel 1991), and an electron temperature 
(Te) of 30,000 K for the CSM. At a SN shock velocity of 
9, 250 km s^^, by 4,300 days after the explosion the front 
has reached a distance in the pre-SN wind, with assumed 
velocity 10 km s~^, where the material was lost from the 
star ^ 11,000 years before the explosion. At such a mass- 
loss rate (M) and wind velocity (w) the CSM, through 
which the shock has already passed, amounts to ^ 1.7M0, 
a significant amount of material, even from a massive pro- 
genitor star. The additional matter present in the next 
layers, i.e. between 4,300 and 7,100 days, corresponds 
to another ~ 2.2Mq, making a total mass ~ 3.9 Mq 
of swept-up CSM by 1998 October. If the flux density 
continues to evolve in the present manner, the amount 
of engulfed matter increases rapidly: for example, in an- 
other 8 years (i.e., by the year 2007) the additional swept 
mass would be ~ 2.9Mq, giving a total swept-up mass of 
~ 6.8 Mq. 

Such a swept wind mass is large, even for a red super- 
giant star, and it implies a very massive progenitor. Es- 
timates of the mass of the envelope ejected by SN 1979C 
range from Menv ~ IM© (Chugai 1985) to M^nv ~ QMq 
(Branch et al. 1981; Bartunov & Blinnikov 1992; Blinnikov 
& Bartimov 1993). In particular, in the models by Bar- 
tunov & Bhnnikov (1992) for SN 1979C, the CSM has a 
reasonable density profile, matching the observed B light 
curves of SN 1979C quite well. Adopting their value for 
Menv, assuming a stellar remnant mass of Mrem ~ 1.4Mq, 
and using our above estimate for the current value of the 
swept-up mass, the initial mass of the progenitor must 
have been Mg ^ 11.3M0. Assuming the current evolution 
of the flux density continues, by 2007 our estimate would 
rise to Mq ^ 14.2Mq. While large, this is still consistent 
with the estimates of Van Dyk et al. (1999). From their 
HST imaging of the SN 1979C environment, they find that 
the stellar ages in that environment are consistent with the 
SN progenitor having an initial mass of Mq w 17-18M0. 

5.2.2. Equatorial wind or disc 

Alternatively, the increase in flux density could be an 
effect of the geometry of the CSM. It is conceivable that 
at small radii the CSM was distributed in a disk with con- 
stant solid angle. At a radius of ~ 3.4 x 10^^ cm (0.11 
pc, the radius reached by a 9,250 km s~^ shock after 4,300 
days) the disk thickness increases so that it covers a larger 
solid angle, reaching about twice the initial solid angle by 
day 7,100 while maintaining a p oc behavior. Such 
a flared geometry would ease the mass requirement if the 
original disk subtended a relatively small angle (O) and is 
at a large inclination angle, so as to be seen almost edge- 
on. In such a case the mass-loss rate determined by radio 
light curve fitting (in which the rate is derived essentially 
from the free-free absorption of the CSM) would only be 
valid within the solid angle subtended by the disk, rather 
than over the full Att steradians. Then, the mass engulfed 
by day 4,300 would be reduced by a factor ft /An and could 
be as small as ~ 1/10 of the estimates in the previous sec- 
tion, without making the requirement on the viewing incli- 
nation angle too severe. However, even in such a case, the 
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flux density flattening observed since day 4.300 requires 
an increase of the interacting surface, so that the mass en- 
gulfed at later times increases at the same rate as the flux 
density does relative to the best-fit model. For example, if 
we take an initial solid angle of the disk to be = 47r/3, 
the mass swept by day 4,300 would be ^ 0.57 Mq and 
the mass swept-up by day 7,100 would be an additional 
~ 0.54 -Mq. More generally, if we denote with (p the frac- 
tion of solid angle subtended by the disk in the inner re- 
gion, the mass swept by day 4,300 would be 1.710 Mq 
and the mass subsequently swept up by day 7,100 would 
be ~ lAAipMQ, a large, but much smaller amount than 
for the spherical shell case. 

5.2.3. Model discrimination 

A possible discriminant between the two basic possibil- 
ities, an increase of CSM density or an increase of CSM 
coverage, is the measurement of free- free absorption at low 
frequencies. For the case of flattening flux density due to 
a less rapidly declining (as p oc r~^'^ rather than p oc r~^) 
CSM density the emission measure will begin to decrease 
as EM cx r~^'^, while for the case of flattening flux density 
due to increasing coverage factor, the emission measure 
behavior will be either the same as, or steeper than the 
canonical stellar wind p oc and could vary as much as 
EM oc or more. In temporal terms, the differing den- 
sity dependence of the EM implies that up to 4,300 days 
the emission measure decreased as EMcx i""^, and that for 
t — to > 4,300 days its behavior changed to EMoc t~^'^. 
However, this behavior can only be tested at those fre- 
quencies where the optical depth was of the order of unity 
at 4,300 days, so that observations at lower frequencies 
than the 1.4 GHz will be required to distinguish between 
the two scenarios. 

To estimate how low an observing frequency is required 
to see the difference, the best-fit model gives an optical 
depth of 

so that at 4,300 days the free-free optical depth was of the 

order of unity at 148 MHz and ~ 0.3 at 250 MHz. Even 
330 MHz observations at the VLA are insufficient as they 
would probe an optical depth of only ^ 0.2 at day 4,300. 
Thus, testing these models may be impractical, since no 
high resolution, high sensitivity radio telescopes currently 
exist at such low frequencies. 

5.3. The Possibility of a Clumpy CSM 

We have previously postulated that the CSM around 
SN 1979C is highly structured (Weiler et al. 1991, 1992), 
and we and others have found evidence for a clumpy CSM 
in SN 1986J (Weiler et al. 1990) and SN 1988Z (Van Dyk 
et al. 1993; Chugai & Danziger 1994). We may also spec- 
ulate that there is evidence for a dense clumpy medium 
surrounding SN 1979C from the fact that the 1.4 GHz 
fiux densities at t > 4,300 days remain relatively constant 
and consistently higher than predictions from the best-fit 
model, while the 5 GHz fiux density values fiuctuate from 
as low as the best-fit model extrapolation, to as high as 



0.47 times the 1.4 GHz fiux densities, with a range of spec- 
tral indices from about —1 ^ a^^ ^ —0.6 (as shown in Fig- 
ure 3). These fiuctuations appear real in that they greatly 
exceed the estimated measurement errors and seem to have 
a time scale of < 1 year. Unfortunately, the paucity of ob- 
servations in the 4,300 7,100 day interval makes it hard to 
test this hypothesis in detail. 

If real, this 5 GHz fluctuation and 1.4 GHz stability 
could indicate a "cooling" time of ~ 0.5 year for the rela- 
tivistic electrons responsible for the higher frequency emis- 
sion and an appreciably greater time constant for lower en- 
ergy electrons. This could be evidence for the presence of a 
relatively small number of dense clumps interspersed in the 
general, stellar wind-generated CSM, but requires a cool- 
ing effect which is a very strong function of frequency. Syn- 
chrotron and inverse Compton losses scale only as 
and Coulomb losses as v^^^ (see, for example, the analy- 
sis of SN 1993J by Fransson & Bjornsson 1998), so that 
producing such a large variation in time scales over such 
a small frequency range by these mechanisms is very dif- 
ficult. Additionally, in order to produce such losses at all, 
synchrotron cooling would require magnetic fields i? ~ 1 
G, much greater than the expected 0(1 mG) fields. Un- 
fortunately, there are too few measurements at 8.4 and 
14.96 GHz (only two at each frequency) to determine if 
the fiuctuations are also seen at higher frequencies. 

Other evidence for dense clumps comes from the esti- 
mate of high densities (rig > 10® cm^'^) derived by Fesen et 
al. (1999) from optical and UV spectroscopy of SN 1979C. 
However, it is hard to assess the significance of their result, 
because their density estimate is based on a comparison of 
[O II] and [O HI] line intensities measured at two epochs 
about 4 years apart, with the unverified assumption that 
the line fitixes do not change with time. Additionally, the 
optical emission and radio emission probably arise from 
different physical media, making it difficult to directly 
compare the densities derived by Fesen et al. (1999) with 
the radio results. Chugai & Danziger (1994) discriminated 
between equatorial disc and clump models for the case of 
SN 1988Z by using the intermediate velocity component of 
the emission lines with FWHM of ^ 2,000 km s^^. Fesen 
et al. (1999) propose that the "spiky" profiles of several 
lines from SN 1979C are suggestive of clumpy emission re- 
gions, and individual spikes in the profiles of [O I]A6300 
and [O II]A7325 in their Figure 5 have roughly the same 
velocity width as the lines used by Chugai & Danziger 
(1994)' for the case of SN 1988Z. 

The possible presence of clumps in the CSM and the 
reality of high frequency fiuctuations should be tested. 
More frequent, multi-frequency radio observations (~ ev- 
ery three months) at a number of frequencies > 1.4 GHz 
might establish the nature of the apparent radio fiux den- 
sity fiuctuations, and simultaneous optical/UV observa- 
tions would permit an unambiguous determination of the 
gas density in the clumps. We have already begun more 
frequent monitoring of SN 1979C with the VLA includ- 
ing both higher (8.4 & 14.9 GHz) and lower (330 MHz) 
frequencies. 

6. CONCLUSIONS 

Analysis of the radio emission from the Type II 
RSN 1979C at 20 and 6 cm from 1991 July 28 through 
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1998 October 19, has shown that its radio emission has 
unexpectedly stopped decreasing in flux density and has 
flattened, or perhaps begun increasing, while maintain- 
ing a relatively constant spectral index. Such behavior is 
in conflict with the best-fit model parameter predictions, 
based on an assumed p oc CSM estabUshed by a con- 
stant mass-loss rate, constant velocity wind from the pre- 
SN star. We interpret this "flattening" to indicate that 
the SN shock wave has encountered a new region of CSM 
which was formed by the SN's progenitor ^ 10,000-15,000 
years before the SN explosion. 

Interpretation of the data implies that this new region 
could either be a higher density shell, which should soon 
be crossed by the fast moving shock, or a "flared" disk-like 
structure in the CSM, if the mass-loss were constrained to 
a narrow solid angle. Additionally, rapid radio flux den- 
sity fluctuations at 5 GHz, which are not present at 1.4 
GHz, are interpreted as possible evidence for clumps or 
large scale density enhancements in the CSM. 

Continued monitoring of SN 1979C at multiple radio 
frequencies - which is ongoing- is needed to determine the 



form and duration of this new phase in the evolution of 
the radio light curves and, correspondingly, the structure 
of this new component of the CSM. With the longest, rela- 
tively complete, multi-frequency data set available for the 
emission from any supernova, SN 1979C serves as a unique 
laboratory for understanding of the evolution of Type H 
supernova progenitors, their pre-SN mass-loss history, and 
their interactions with their local environments. 

We wish to thank the referee, Stephen Reynolds, for 
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MJM wish to thank the Office of Naval Research (ONR) 
for the 6.1 funding supporting this research. CKL addi- 
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and Engineering Apprenticeship (SEAP) program at the 
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Table 1 

Measured flux density values for the secondary calibrator 1252+119 



Time Since 



Observation 


SN 1979C 


S20 


^6 

(Jy) 




^2 


Date 


Optical Max. (days) 


(Jy) 


(Jy) 


(Jy) 


1979 Apr 19 


= 










1991 Jul 28 


4483 


0.770 


0.607 






1991 Oct 31 


4578 


0.834 


0.717 






1992 Mar 01 


4700 


0.776 


0.628 






1992 Oct 13 


4926 


0.757 


0.674 






1993 Jan 28 


5033 


0.795 


0.672 






1993 May 07 


5132 




0.673 






1993 Oct 17 


5295 


0.771 


0.717 






1994 Feb 18 


5419 


0.789 


0.738 






1994 Apr 25 


5485 


0.753 


0.721 






1995 Jun 15 


5901 


0.752 


0.677 






1995 Dec 12 


6081 


0.771 


0.746 






1996 Oct 06 


6380 


0.772 


0.772 






1996 Dec 21 


6456 


0.736 


0.791 


0.753 


0.704 


1997 Sep 23 


6732 


0.740 


0.702 






1998 Feb 10 


6872 






0.845 


0.876 


1998 Feb 13 


6875 


0.769 


0.720 






1998 Oct 19 


7123 


0.785 


0.771 
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Table 2 

New Flux Density Measurements for SN IQTQC* 



Obs. 
Date 



Time 

Since 
Optical 
Max.'' 
(days) 



Flux Density (mjy) 



VLA 
Conf. 



S20 



C20 



(T6 



S3.6 C3.( 



(T2 



1979 


Apr 19 


= 




















1991 


Jul 28 


4483 


A 


5 


229 





276 


1 


915 





103 


1991 


Oct 31 


4578 


B 


5 


365 





289 


1 


546 





101 


1992 


Mar 01 


4700 


C 


5 


640 





801 


2 


253 





120 


1992 


Oct 13 


4926 


A 


5 


396 





287 


2 


574 





143 


1993 


Jan 28 


5033 


A 


5 


582 





304 


2 


625 





136 


1993 


May 07 


5132 


B 










2 


583 





138 


1993 


Oct 17 


5295 


C/D 


5 


190 





281 


2 


251 





115 


1994 


Feb 18 


5419 


D/A 


5 


670 





535 


1 


660 





311 


1994 


Apr 25 


5485 


A 


5 


810 





314 


2 


740 





170 


1995 


Jun 15 


5901 


A 


5 


340 





328 


1 


740 





240 


1995 


Dec 12 


6081 


B 


5 


663 





074 


2 


664 





139 


1996 


Oct 06 


6380 


D/A 


5 


640 





335 


2 


170 





170 


1996 


Dec 21 


6456 


A 


5 


820 





315 


2 


690 





180 


1997 


Sep 23 


6732 


CnD 


4 


200 





452 


1 


727 





107 


1998 


Feb 10 


6872 


D 


















1998 


Feb 13 


6875 


DnA 


5 


460 





279 


2 


680 





149 


1998 


Oct 19 


7123 


B 


5 


140 





406 


2 


500 





136 



:.ooo 0. 

.930 0. 



112 
107 



1.120 
0.950 



0.195 
0.264 



''For previous measurements, cf. Weiler et al. (1986, 1991). 

''The date of the explosion is taken to be 1979 April 4, 15 days before optical maximum (cf Weiler 
et al. 1986). 



Montes et al. 



9 



Table 3 

Fitting Parameters for SN 1979C 



Parameter 


Value 


Deviation Range^ 


Xi(mJy) 


1710 


1440-2060 


a 


-0.75 


-(0.76-0.63) 


P 


-0.80 


-(0.83-0.78) 


K2 


3.38 X 10^ 


(2.97-3.82)xl0^ 


S = a-f3-3 


-2.94 


-(2.96-2.92) 
(6.7-8.1)xl0-2 


A 


7.3 X 10-2 


B-i(days) 


1570 


1520-1610 


C 


O.gOTT 


(0.83-0.99)7r 


to 


=1979 Apr. 4 




xVdof 


2.47 





'^Only data through 1990 December are used in the 
fitting procedure. 

'^The error estimates for the parameter values are 
determined using the bootstrap method, which is de- 
scribed in §11 
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1252+1 19 




1000 



2000 3000 4000 5000 

Days since 1979 April 4 



6000 7000 



Fig. 1. — Measured flux densities for the secondary calibrator 1252+119 at 6 cm {filled triangles) and 20 cm {open circles). 
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V. 




500 1000 2000 

(t-to)/days [t^^ (19-Apr-79) -15.000 days.] 



5000 



Fig. 2.— Radio "light curves" for SN 1979C in NGC 4321 (MlOO) at wavelengths 20 cm (crosses) and 6 cm (open 
circles). (Since there are only two new measurements available at 3.6 and 2 cm, they are not plotted here. However, 
they are shown as confirming the nonthermal nature of the spectrum in Figure ^.) The data represents about 18 years of 
observations for this object, including the new observations presented in this paper and previous observations from Weiler 
et al. (1986, 1991). The curves represent the best-fit model light curves at 6 cm (solid), and 20 cm (dashed), including 
the quasi-periodic, or sinusoidal, term proposed by Weiler et al. (1992). The best-fit parameters were determined using 
only data through 1990 December (day ^ 4,300). 
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^ — I — I — I — — I — I — \ — I — — I — I — I — I — — I — I — I — I — — I — I — I — I — — I — h 




1000 2000 3000 4000 

(t-g/days [to= (19-Apr-79) 



5000 6000 
15.000 days.] 



7000 



Fig. 3. — Evolution of spectral index a (where S oc between 20 and 6 cm for SN 1979C, plotted as a function of 

time (in days) since the estimated explosion date of 1979 April 4 (15 days before optical maximum). The solid line is 
calculated from the best-fit theoretical "light curves" shown in Fig. ^ For reference, see Figure 3 in Weiler et al. (1986) 
and Figure 2 in Weiler et al. (1991). 
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Fig. 4. — The radio spectrum for SN 1979C on 1996 December 21 {open circles) and 1998 February 10 & 13 {open stars) 
from 20 to 2 cm. Note that, to within the errors, the spectra are completely optically thin at all frequencies. The solid 
line is the best-fit spectrum to the combined data from both dates. The best-fit spectral index value for the combined 
data on these dates is a = —0.63 ± 0.03. 



